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Abstract 

We  present  a  mathematical  framework  to  provide  general  well-posedness  and 
approximation  results  for  a  laminated  curved  beam.  The  layered  beam  consists 
of  an  elastic  core  layer,  symmentric  viscoelastic  damping  layers,  symmetric 
constraining  layers  and  bonded  piezoceramic  patch  pairs  on  the  outer  surfaces. 
The  model  includes  unbounded  inputs  from  the  piezoceramic  patches,  hysteresis 
and  shear  in  the  viscoelastic  layers.  The  viscoelastic  layers  are  held  in  place 
by  sandwiching  them  between  the  elastic  core  and  elastic  constraining  layers. 
Families  of  linear  and  cubic  splines  are  used  to  illustrate  the  approximation 
ideas. 
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1  Introduction 


There  is  a  widespread  interest  in  the  engineering  and  industrial  community  over 
the  past  several  decades  in  the  use  of  composite  (elastic/viscoelastic  smart  mate¬ 
rial)  structures  such  as  beams,  plates  and  complex  articulated  structures.  There  is  a 
large  body  of  literature  dealing  with  the  design  of  viscoelastic  material  (VEM)  lay¬ 
ers  to  enhance  passive  damping  ([MM],  [YD],  [R],  [.JK] ,  [Ma],  and  references  therein). 
Recently,  however,  there  is  a  growing  interest  in  active  constrained  layer  (ACL)  struc¬ 
tures  ([ATW],  [HWT] ,  [RWTW],  [RWTW1],  [LW],  [KKS],  [BGM],  [BMZ],  [BR],  [VKI] 
and  references  therein).  In  such  structures  sensor /actuator  devices  such  as  piezoce¬ 
ramic  patches  or  layers  with  appropriate  circuitry  and  power  supply  are  bonded  to 
a  parent  material  that  is  elastic.  Between  the  parent  material  and  the  piezoceramic 
patches  another  material,  VEM  is  often  embedded  to  provide  passive  damping. 

There  is  very  little  literature  on  rigorous  mathematical  formulations  for  the  ACL 
devices  mentioned  above.  The  present  paper  is  a  continuation  of  the  simpler  model 
in  [BMZ],  Here  we  have  added  constraining  layers  to  hold  the  viscoelastic  layers  in 
place,  and  further  we  account  for  the  shear  in  the  viscoelastic  layers.  As  in  [BMZ] 
we  establish  well-posedness  and  numerical  approximation  framework.  The  multilayer 
curved  beam  configuration  presented  produces  fully  coupled  transverse  and  longitu¬ 
dinal  vibrations.  In  addition  to  the  shear  we  have  included  “integral  hysteresis”  in 
our  VEM  layers.  Such  curved  configurations  with  VEM  hysteresis  are  important  in 
automobile  and  aircraft  industrial  applications. 

Our  approach,  which  is  functional  analytic  in  nature,  combines  the  ideas  on  treat¬ 
ment  of  time  hysteresis  and  approximation  in  partial  differential  equation  models 
given  in  [BFW1]  and  [BFW2]  with  general  well-posedness  (in  terms  of  semigroup) 
arguments  given  in  [BIW]  and  [BKW].  These  ideas  and  results  can  be  generalized 
and  extended  to  treat  many  of  the  more  realistic  (linear)  physical  models  found  in 
the  literature.  The  resulting  theoretical  and  computational  framework  can  be  used 
in  control,  estimation  and  simulation  applications. 

In  the  next  section  we  present  our  example  model  along  with  a  brief  outline  of  the 
elements  of  its  derivatives.  We  recast  this  model  in  abstract  Cauchy  system  form  in 
Section  3  and  show  that  it  can  be  associated  with  a  Co-semigroup  in  appropriately 
chosen  spaces.  In  Section  4,  we  give  convergence  results  for  a  typical  family  (linear 
splines  for  longitudinal  vibrations,  cubic  splines  for  transverse  vibrations)  of  approxi¬ 
mations.  We  have  successfully  used  this  family  in  related  computational  studies,  the 
results  of  which  are  detailed  elsewhere  (see  [BZ]). 


2  Model  Development 

2.1  Physical  Model  and  Assumptions 


We  study  a  laminated  beam  that  contains  an  elastic  core  with  two  viscoelastic 
layers,  two  constraining  layers,  and  two  piezoceramic  patches.  The  viscoelastic  layers 
are  bonded  to  the  elastic  core,  the  constraining  layers  to  the  viscoelastic  layers,  and 
the  piezoceramic  patches  to  the  constraining  layers.  Pairs  of  the  viscoelastic  layers, 
constraining  layers,  and  piezoceramic  patches  are  used  to  produce  geometrical  sym¬ 
metry.  Each  component  has  a  finite  radius  of  curvature.  A  view  of  a  middle  section 
of  the  composite  curved  beam  from  the  side  is  given  in  Figure  1. 

The  transverse  vibrations  of  the  composite  beam  are  our  primary  consideration. 
However,  due  to  the  beam  curvature  there  are  two  displacements  involved:  v{t,6J  in 
the  tangential  (6)  direction  and  w(t:6)  in  the  radial  (h)  direction.  Any  motion  in 
the  direction  of  the  width  of  the  beam,  x,  is  assumed  negligible.  We  assume  that 
the  radial  displacement  w(t,  0)  is  the  same  for  all  layers  and  we  derive  tangential 
displacements  for  each  layer.  Moreover,  we  make  the  following  assumptions  for  our 
model: 

•  The  viscoelastic  material  under  consideration  is  linear. 

•  Every  material  layer  is  bonded  perfectly  to  others  in  contact,  i.e. ,  no  slip,  and 
no  disbonds  are  allowed. 

•  Every  material  layer  deforms  uniformly  with  respect  to  thickness;  that  is,  dis¬ 
placements  are  only  functions  of  t  and  6. 

•  The  thickness  of  the  composite  beam  is  very  small  in  comparison  with  other 
dimensions  such  as  radius  of  curvature  and  length. 

•  Composite  beam  deformations  are  sufficiently  small  so  as  to  allow  higher  order 
motions  to  be  neglected,  i.e.,  we  use  first  order  models. 

•  Transverse  normal  stresses  (i.e.,  the  stresses  in  the  direction,  h,  normal  to  the 
thin  dimension)  are  taken  to  be  negligible. 

•  A  cross  section  which  is  originally  normal  to  the  reference  surfaces  of  the  elastic 
core,  the  constraining  layers,  and  the  piezoceramic  patches  will  remain  normal 
to  the  deformed  reference  surfaces  and  will  remain  unstrained. 


E: Elastic  Core  thickness  h 
VE:  Viscoelastic  layer,  thickness  A 
C :  Constraining  layer,  thickness  r 
P :  Piezo  patch,  thickness  T. 

Figure  1.  Middle  Section  of  the  Curved  Composite  Beam 

2.2  Mathematical  Development  of  Model 

To  develop  the  mathematical  model  we  treat  each  of  the  layers  as  a  shell.  Fur¬ 
ther,  we  introduce  four  sublayers,  h8j,  i  —  1,2,  j  —  1,2,  with  nonzero  thickness  A *j, 
respectively  in  the  viscoelastic  layers  (see  Figure  2).  We  would  like  the  viscoelastic 
layers  to  have  the  same  circumferential  displacements  as  the  elastic  core  at  the  inter¬ 
faces  between  the  two  materials.  We  would  like  a  similar  equality  at  the  interfaces 
of  the  viscoelastic  layers  with  the  constraining  layers.  However,  unlike  the  elastic 
and  constraining  layers,  the  viscoelastic  layers  undergo  shear  distortion.  This  can 
be  efficiently  captured  by  considering  sublayers  of  the  viscoelastic  layer  divided  by  a 
neutral  line  of  zero  shear  which  is  no  longer  necessarily  the  geometrical  center  line. 
Moreover,  its  location,  which  changes  with  time  during  deformations,  is  in  general 
unknown.  By  introducing  equal  ions  for  two  sublayers  in  each  viscoelastic  layer,  we 
can  account  for  this  neutral  shear  axis  without  actually  knowing  its  location.  We 
emphasize  that  this  is  simply  a  device  to  develop  a  model  which  allows  shear  in  the 


viscoelastic  layers  as  well  as  different  materials  in  the  constraining  layers  and  elastic 
core. 


The  viscoelastic  layers  divided  into  sublayers  in,  i12,  T21  and  L2 1- 
Figure  2. Viscoelastic  Layer  Subdivisions 


In  what  follows  we  write  dt  for  d/di  and  d  for  d/d6.  In  the  elastic  core,  referring 
the  motion  to  its  neutral  line  we  write 


>(t,  9 ,  z)  =  v(t,  9)  +  z 


'v(t,9)  dw(t,9) 


R 


R 


h  h  , 

"2  -  Z  £  2'  (2'1) 


where  R  is  the  radius  of  curvature,  h  the  thickness.  The  formula  (2.1)  comes  from 
standard  shell  theory  (see  [BSW]).  The  strain  at  any  point  on  the  elastic  core  is  given 

by 

ee(t)  =  °)  +  °)  ~  ^2  °)  (2-2) 

and  the  corresponding  stress  is  given  by 

°e{t)  =  Y^ee{t),  (2.3) 


where  Ee  is  the  Young’s  modulus,  and  ve  the  Poisson  ratio.  In  (2.1)  and  (2.2)  the 
variable  z  is  the  normal  distance  of  the  point  (61,  z)  from  the  neutral  line. 

Moving  up  to  the  viscoelastic  layers  In,  and  Ly2  we  treat  each  sublayer  as  a  shell 
and  refer  its  motion  to  its  geometrical  center  line.  Since  we  include  shear  deformation 
in  these  layers  (Tjj),  Kirchhoff’s  fourth  hypothesis,  which  implies  7^  =  7^  =  ez  =  0 
is  no  longer  valid.  Denoting  the  circumferential  displacements  by  VM(t,  9,  z),  and 
(/.,  9 ,  z)  we  write 
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11 


<  5  < 


A 


11 


v1v?(tA*)  =  v£{t,e,o)  +  z8n, 


2 


2 


(2.4) 


fefeMfe  =  v™(t,6,0)  -zSJ2,  <  5  <  (2.5) 

In  (2.4)  and  (2.5),  5  denotes  the  normal  distance  of  a  point  ( 6 ,  z)  from  the  geometrical 
centerline.  The  fact  that  Kirchhoff’s  fourth  hypothesis  is  not  valid  in  the  sublayers 
Lij  prohibits  us  from  expressing  StJ  in  terms  of  Lfe(f,  6,  0),  w,  and  RtJ.  the  radius  of 
curvature  of  the  geometrical  center  line  of  the  sublayer  Lij ,  in  a  simple  relationship 
of  the  form  (2.f).  Thus,  we  have  to  seek  an  alternative  way  of  determining  the  StJ 's. 
This  will  be  done  later. 

To  express  the  strain  fe'fe/.,  Of)  and  the  stress  (aget j(t ,  Of)  relationship  in  the 
sublayer  Lij  we  use  Boltzmann  superposition  model.  Thus, 

W  =  T—T  fejW  -  f  sW%(‘  +  «)*)  .  (2-6) 

1  -  V  j-r  J 

where 

<%lt)  ~  2-  (dV» (l,  e ,  z)  +  w(t , «))  ,  (2.7) 

and  Eve  is  a  modulus  of  elasticity,  vve  the  Poisson  ratio  of  the  material. 

To  express  the  shear  strain  (  7<;1  feh  Of)  and  shear  stress  0)  relationship  we 

again  use  the  Boltzmann  superposition  model.  Thus, 

fefefe)  =  2(1  +  ^  pfefe)  -  J_r  +  s)ds\  «  (2-8) 

where 

=  irfe  («»«.*)  -  +  (-l)‘+)S.r  (2.9) 

Next,  we  move  up  to  the  constraining  layer.  Let  Ufefez)  denote  the  circum¬ 
ferential  displacement,  and  the  radius  of  curvature  of  the  neutral  line.  Then,  we 
have 


Ife/feife-  fe(/fe,0)  •  5 


fe(/fe,0)  chr(/.,fl) 


(2.10) 


where  r  is  the  thickness  of  the  constraining  layer  and  z  is  the  local  coordinate  for  the 
distance  from  its  neutral  line.  Using  a  form  similar  to  (2.2)  for  the  elastic  core  layer, 
the  strain  is  given  by 

7i(‘. «)  =  f  (dvnt,  0-  0)  +  »(i, «))  -  «)■  (2-H) 

The  strain  (egfet,  #))  and  stress  (cfefe#))  relationship  is  given  by 

<!(«)  =  (2-i2) 


where  Ec  is  the  Young’s  modulus,  ue  the  Poisson  ratio. 


Finally,  we  move  up  to  the  piezoceramic  patch.  Let  Vf(t,  9 ,  z)  be  the  circumfer¬ 
ential  displacement,  the  radius  of  curvature  of  the  neutral  line,  T  the  thickness. 
Then,  we  write 


Vjp(t,0,z) 


'WM,o) 


dw{t ,  0)  \ 

r )i 


T  T 

- <  5  <  -. 

2  “  “  2 


(2.13) 


Once  again  the  strain  is  given  by 


«!,(*)  =  ^,*0)  +  ^-^ 


d2w. 


The  strain  (e^-^i))  and  stress  (a:0  ^  (I.))  relationship  is  given  by 


Ep 


<1  (0  = 


(2.14) 


(2.15) 


where  Ep  is  the  Young’s  modulus,  vp  the  Poisson  ratio,  and  Xi(@)  =  1  where  the  piezo 
patch  is  present  and  zero  elsewhere. 

We  repeat  the  above  steps  for  the  layers  below  the  elastic  core.  Starting  from  the 
elastic  core  we  move  down  to  encounter  the  sublayer  L2 1,  then  T22.  Let  V££ ,  i  =  1,2 
be  the  circumferential  displacements  of  these  sublayers  relative  to  the  geometrical 
center  lines.  Then,  we  write 


V£{t,9,z)  =  V£{t,9,Q)-~zS21, 


v-(M,5)  =  y2T(M,0)  +  5,s22, 


A 


21 


2 

A22 


^  ~ .  -  ^21 

<  Z  <  - , 

-  -  2 

(2,16) 

a22 

<  z  <  - . 

~  “  2 

(2.17) 

Moving  down  to  the  constraining  layer,  we  let  V2[i.  9 ,  z )  denote  the  circumferential 
displacement,  and  lir2  the  radius  of  curvature  of  the  neutral  line.  Thus,  we  write 


v;{t..9.i) 


v2c{t,9, 0)+  ! 


(VfjtA  0) 


dw(t,  0)\ 


T  T 

—  <  5  <  — , 
2  -  “  2’ 


(2.18) 


where  r  is  the  thickness  of  the  constraining  layer.  This  constraining  layer  interfaces 
the  viscoelastic  layer  above  it  when  5  =  r/2  and  the  piezo  patch  below  it  when 
z  =  —  t/2. 

Finally,  we  move  down  to  the  piezo  patch.  Let  V2 ( i .  9 ,  z)  be  the  circumferential 
displacement,  ll2  the  radius  of  curvature  of  the  neutral  line,  T  the  thickness.  Then, 
we  write 


Vj (t,  9 ,  0)  +  5 


fgM  Q) 

V  Rl 


dw{t ,  #)\ 

m  ) 


T  T 

- <  5  <  -. 

2  “  “2 


(2.19) 


When  5  =  T/2  the  piezo  patch  interfaces  the  constraining  layer  above  it. 


Again,  as  we  did  for  the  layers  above  the  elastic  core,  we  write  the  stress-strain 
relationships.  In  the  viscoelastic  layers  we  have  shear  displacements  and  have  ignored 
compression. 

Finally  the  no-slippage  and  no-disbond  conditions  allow  us  to  give  the  circumfer¬ 
ential  displacements  and  the  Su*p  (see  (2.4),  (2.5),  (2.16),  (2.17))  in  terms  of  v:w, 
and  h)je(^,^,0).  Thus,  we  take  v ,  w,  and  V) )e(f,  6,  0)  as  our  dynamic  state  variables. 


2.2.1  Internal  Forces,  Moment  Resultants  and  Equations  of  Motions 


In  the  elastic  layer,  constraining  layers,  and  piezo  patches  we  denote  the  total 
internal  force  by  tV|,  Ng^  i  =  1,2,  and  Ng  i  =  1,2  respectively.  Here  Ngt  and 
Nh  are  the  total  internal  force  resultants  along  the  neutral  surfaces  in  the  upper 
constraining  layer  and  upper  piezo  patch.  Similarly  Njj2  and  Ng2  denote  the  total 
internal  force  resultants  for  the  lower  ones.  Next,  the  total  internal  force  resultant 
along  the  geometrical  mid  surface  in  the  viscoelastic  layer  Li j  is  denoted  by  Ay 
i  =  1,2,  j  =  1,2.  The  corresponding  resultant  shear  forces  in  the  viscoelastic  layers 
are  denoted  by  QYtJ-  i  =  1,2;  j  =  1,2. 

We  have  (see  [BSW],[BSW2] ) 

m  = 


where 


and 


<T, 


m  = 


Continuing,  we  have 

r!2 


Nh  = 


/ 
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(To  idz 

-r/2  ’ 

/2  ec 


rh/2 

/  ° 

J-h/2 

ei  Eeh 

'ldz=  fi(l > 

(2.20) 

Nve 

/■Asj/2 

(2.21) 

Eve  1 

<  o.ri1)  -  j  9{s)eVe%{t  +  s)ds 

J  —r 

,  (2.22) 

1 

0) 

*  4-  w + »)■ 

(2.23) 

(2.24) 

-r/2  1  - 


^dW(tA0)  +  ^w--(^f^w)j  dz , 


NL  = 
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P  :d,Z 


-T/2 

T/2  Ep 


-T/2  1  —  vl  \R 


—  dv?(t,e,o)  +  —w- 
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(R: 


•P\2 


d2w 


The  shear  in  the  viscoelastic  layer  is  given  by 


Qe%  = 


L 


A.j/2 


-A,j/2 


x!m 


z. 


(2.25) 


(2.26) 


where 


<h+)  = 


E’ 


2(1+  I+e) 


7  «.,,(<)  -  1°  g(*h, ;;,«(« +  3)* 

«/—  r 


Next  we  present  the  internal  moments 
A/2 


Ml  = 
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(Jgzdz 

—  h/2 

Ee  rh/2  ( dv 


(2.27) 

(2.28) 

(2.29) 


iA  t  dv  w  z  r 

1  -  1/2  i_ft/2  +  E  _  h*2  d 


w  z 


~Ee}l3  02 
«z  =  - —  d  w. 


l-h/2 


Y2R2{1  -  v2 ) 

The  internal  moments  for  the  constraining  layers  about  the  neutral  line  of  the 
elastic  core  are  given  by 

Ec  r>2  (  1 


ME  = 
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1  -  v2  J-r/ 2  \R 

tEc  W  -  R 


dV?(t,  0,0)  +  —  -  — —  d2w  (R/  +  5  -  R) dz 


m  m 


1  -  m2  Ri 


{dVf(t,  0, 0)  +  w)  - 


SRC  Q2 


W 


12(1  —  i/2)  {Rfy 


Thus, 


MSj  +  MS,  = 


Si 


u=j 
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tEc  Ri  -  R 


—  V2  RC: 


(dv;(L0,0)  -tr) 
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■  d2w 
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(2.30) 


Note  that 


12(1  —  v2)  \(Ri)2  (R$ 

Ri  —  R  ^  h/2  +  An  +  A 1 2  +  t f 2 
Ri  ~  R  ’ 

i?2  —  R  ~{h/2  +  A21  +  A22  +  t /2) 


m 
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The  internal  moments  for  the  piezo  layers  about  the  neutral  line  of  the  elastic  core 
are  given  by  (see  [BSW],[BSW2]) 


ML  = 


tEp  R/-R 

1  -  M 


( dVp(t,9,0)  +  w )  - 


2EPT3  Xm 

3(1  -wm2 


d2w. 


Thus, 
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_  2EpT 3  2,,,/xg  . 

3(i  -vi)  \m2  m2)' 

To  simplify  our  presentation  we  define  the  following  entities: 
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0  a22; 


Die)  =  T{Jve  +  T2{Jc  +  .P) 


(2.34) 

(2.35) 

(2.36) 

(2.37) 

(2.38) 

(2.39) 

(2.40) 

(2.41) 

(2.42) 

(2.43) 

(2.44) 
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3(1 -  vi)  \m2  {my 


(2.45) 


(2.46) 


+ 


(2.47) 


S  —  (Nil,  *5i2,  S2I1  S22)T ■ 

Furthermore,  define  the  resultant  moment  for  the  nonviscoelastic  materials 

Me  =  Ml  +  £  MSti  +  E  Ml .  (2.48) 

i  i 

Kirchhoff’s  fourth  hypothesis  implies  7 xz  =  7 gz  =  ez  =  0  and  hence  we  cannot  use  7 ez 
to  compute  shear  stress  or  the  associated  moment.  However,  from  moment  balance 
one  can  relate  internal  moment  to  shear  stress.  Let  Qg  be  the  sum  of  the  shear  stresses 
in  the  elastic,  constraining,  and  piezo  layers.  Then,  we  have  the  relationship 

Qgn^-dMg.  (2.49) 

H 

Let 

0“  =  £  Q7fl  «  «2(S»  -  v)  +  (Du  S).  (2.60) 

ij 

Thus,  the  sum  Qg  of  all  the  shear  stresses  is  given  by 

Qe  =  Qe  +  Q  T  ■ 

Next,  we  assemble  the  internal  forces  of  the  beam.  Thus,  the  sum  of  the  internal 
force  resultants  Ng  is  given  by 

N,  =  £  iV,”  +  £  Ntf  +  £  Nl,  +  NS. 

ij  i  i 

Here  we  are  treating  the  entire  curved  laminated  beam  as  a  shell  and  Ng  as  the  total 
internal  resultant  force  along  the  middle  surface  of  the  elastic  core. 

Let  p  be  the  linear  mass  density  which  varies  with  6  due  to  the  presence  of  the 
patches.  If  we  take  pe ,  pve ,  pc ,  pp  be  the  volume  densities  of  the  elastic,  viscoelastic, 
constraining,  and  piezoceramic  layers  respectively,  then  we  have 

p  =  p'hb  +  pve2Ab  +  pc2rb  +  pQTbxl  +  TbX  f) 

2A  =  £  Ay. 

ij 

The  tangential  and  radial  displacements  v  and  w  of  a  curved  beam  are  given  by 
[BSW],[BSW2],  [M]  by  the  equations 

pdtv  -  4  dNg  -  ^Qg  =  bqg  (2.51) 

ft  ft 

pd^iv  +  4  -  4  dQe  =  bqn  +  4 dm 


(2.52) 


where  b  is  the  width  of  the  beam.  We  shall  supplement  the  system  (2.51)  and  (2.52) 
by  a  set  of  equations  for 

V  =  (V11,V12,V21,V22)t,  (2.53) 

where  (see  (2.4),  (2.5),  (2.16),  (2.17)) 


%(M)  =  yr<(,«,  o) 


(2.54) 


We  consider  the  layer  l,tJ  for  which  the  equation  of  motion  locally  is  given  by 


Rij  ~\~  z 


davg%  +  — , 

Rij  +  ~ 


ve  ,  o  tie 

aesm  "t  OzO g-^j. 


(2.55) 


The  equation  (2.55)  is  obtained  by  taking  an  infinitesimal  element  and  writing  force 
balance  [M],  Again  we  remark  that  5  measures  the  normal  distance  of  the  point  (0,  z ) 
from  the  geometrical  centerline  of  the  layer  LtJ. 

In  (2.55)  set  5  =  0  to  obtain  the  equation 
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2(1  VVe) 

'  dw{t  +  s,  6)  Vij(t  +  s:6 ) 
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%  +  (-1)8+J% 

+  (-i  y+jSij(t+s,0))  ds 


‘*3 


R, 


Let 


A  = 


Ev 


1  —  v 2 


-d/iagi  R^  1  ,  P12  ,  i?2i  ?  ^22  ) 


X  = 


ve 

Eve 


(  p—  2  p— 2  p—  2  p— 2\T 
A-^ll  5  ^12  5  ^21  5  ■rt22  J 


1  —  Z/2 

P  =  <P’ag(l,  —1,  —1,  1). 

Then,  using  (2.57),  (2.58),  (2.59)  we  rewrite  (2.56)  as 

pved2V  =  d2AV  +  dwL  —  f°r  g(s)(d2AV  +  dwL)ds 


(2.56) 

(2.57) 

(2.58) 

(2.59) 


+  “  Vve^ 


dwL  -  AV  +  PAS 
f  g(s){dwL  —  AV  +  P AS)ds 

J  —r 


(2.60) 


In  (2.60)  V  is  as  in  (2.53)  and  S  =  (-Sii,  P12,  *S  1 ,  $22)T • 

Recalling  that  V)j(l,  9)  =  Vje(t,0,  0)  and  using  the  no-slip  condition  we  have 


Let 


Vl2  —  v  +  7777  v  ~  7777  dw  +  Ana'll  —  —  A12.S12 

A  jEL  A  ft  A 

V21  =  v - v  -\ - dw - A21.S21 

2R  2R  2 

V22  =  ^  —  7777  v  +  — —  dw  —  Aji^i  +  7-  A22  S22  • 

A  jft  A  it  A 


~  /  A  h  h  h  ' 

'  1  “  \1  +  AR,1  +  2^,1_AR’1_AR 


J2  =  — 


h  hhh 
2R'~2R'2R'2R 


Then,  we  rewrite  (2.61)  as 


V  =  Ti5  +  r.i|  +  dw.J2. 


Observing  that  T\  is  invertible  due  to  our  basic  assumption  that  A 4-j  p  0, 
have 

s  =  T^V  -  vl\  '  J\  -  dw'l\  \l2. 

Next,  let 

p  =  Pat-1  -  a 
j2  =  l-  Pat-1.),. 

Using  (2.65),  (2.66),  and  (2.67)  we  rewrite  (2.60)  as  follows 

pved2V  =  d2AV  +  dwL  -  f  g{s){d2AV  +  dwL)ds 

J —r 

(1  -  vm)[PV  -  vPAT-Kp  +  dw.j2 


-  f°  g(s)(PV  -  vPA Tf1  Ji  +  dw.J2)ds\. 

J  —r 

From  (2.51)  and  (2.52)  finally  we  have  the  equations 

pd2v  -  ^d[(pdv  +  w)  +  (2d2w  +  ( D(0 ),  dS )] 

+  —  d  f  g(s){cii(dv  +  w)  +  b\d2  w  +  (Jve :TidS)}ds 

it  J —r 

-  ^[a2(dw  -  v)  +  (D1:S)  +  ^dMg\ 

b 

+  j;  J  g(s){a2(dw  -  v)  +  (Di,  S)}ds  =  bqe 


pd2w  +  —  [£l(3u  +  w)  +  (2d2w  +  {D(6),  35)] 

ti 


(2.61) 

(2.62) 

(2.63) 

(2.64) 
we  then 

(2.65) 

(2.66) 

(2.67) 

(2.68) 


(2.69) 


—  —  f  g(s){ai(dv  +  w)  +  bjd2w  +  (Jve,  T1dS)}ds 

it  J —r 

-  ^d[a2(dw  -  v)  +  (D^S)  +  -^dMg] 

H  ti 

+  j;  j  g{s)d{a2(dw  -  v)+ (D1,S)}ds  =  bqn  + -^dmx.  (2.70) 

The  equations  (2.68),  (2.69),  (2.70)  give  a  system  of  equations  governing  the 
motion  of  the  beam  in  terms  of  the  variables  v,  w.  V .  Note  that  in  (2.69)  and  (2.70) 
the  variable  vector  S  is  to  be  replaced  according  to  (2.65). 

We  rewrite  (2.69)  and  (2.70)  by  substituting  for  S.  First  we  let 


d(6)  =  i>n  ’(t;)  'd{o) 

(2.71) 

6  =  bR^Ci  ~  (D{0)Jj 

(2.72) 

a1  =  bR-1(a1-{Jve,j1)) 

(2.73) 
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(2.74) 
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rH 

(2.75) 

=  bR-\t2-a2)  +  {D1- D{6),J2) 

(2.76) 

bR-1({Jv%J2)-(a2  +  b1))  +  {D1J2). 

(2.77) 

Then. 
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pdt 
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-  d^\dv  +  (bR  1a2  +  { /h  ?  Ji))v 
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+  (D(e)-D1,dV) 

+  f°r  g(s){(a  1  +  a2)3u  —  bR~1a1w  +  D3d2w 
+  (ZJj  -  bR^J™,  3V)}ds  =  +  jjdmx. 


(2.78) 


(2.79) 


The  system  we  should  consider  consists  of  the  equations  (2.68),  (2.78),  and  (2.79). 
We  ignore  the  underlined  terms  in  (2.78)  (see  REMARK  3.1). 

The  system  (2.68),  (2.78),  and  (2.79)  can  be  rewritten  in  the  form 
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1  v(l.  ■  x.  0)  \  ^  / 
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\V(t  +  ,,t))l  <> 


(2.80) 


Using  (2.68)  we  consider  the  equation 


pved2tz  =  d2A Z  +  d2wl  -  f  g(s)(d2AZ  +  d2wL)ds 
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-  f  g(s)(PZ  -  dvPAT^-h  +  d2w.J2)ds. 

J —r 

In  (2.78)  and  (2.79)  we  replace  V  according  to  the  formula 

dv(t,0)  =  z(t,e ) 

and  obtain  a  system  corresponding  to  (2.80): 
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(2.81) 


(2.82) 


(2.83) 


where  the  operators  Ai,i7i  (which  are  rather  tedious  to  define)  are  given  in  the  Ap¬ 
pendix. 


Next,  let 


A  —  A\  T  oiB\ 


where 
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a  —  /  g[s)ds. 

J —r 


Then,  we  rewrite  (2.83)  as 
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where  c((.s)  =  v(t  +  s),  wt(s)  =  w(t  +  s),  Zt(s )  =  Z{t  +  s). 


(2.84) 


2.2.2  External  Forces  and  Moments 


We  assume  for  our  model  that  the  external  forces  and  moments  are  of  a  form  such 
as  those  produced  by  piezo  materials.  Thus,  we  ignore  any  contribution  from  external 
air  damping.  The  PZT  patches  can  be  used  as  both  sensors  and  actuators,  and  as 
actuators  they  can  generate  significant  external  forces  and  moments.  We  follow  the 
formulations  in  [BSW],  and  to  aid  in  quantitative  description  of  the  PZT  forces,  we 
define  the  indicator  function 


f  1  9  >  (01  +  e2)/2 

Si,2(0)  =  <  0  6  =  ($1  +  02)/ 2  . 

l-l  0<{01  +  02)/2 


Then,  assuming  that  the  two  piezo  patches  have  the  same  Young’s  modulus,  Ep , 
the  same  Poisson  ratio,  vv ,  and  the  same  strain  constant  d31 .  the  total  extended  force 
from  the  piezo  patches  is  given  by 

Nl  =  (Nhxt  +  Nlzx^ 

=  ~-r~Z — +  V2X2)Si,2- 

1  ISp 

Here  V\  and  Vi  are  the  voltages  applied  to  the  patches  to  produce  deformations. 
Similarly,  the  total  external  moment  due  to  the  patches  is  found  to  be 


=  +  f 


2  +  T  +  A  +  r 


-  +  A  +  r 
2 


Finally  the  external  forces  in  (2.51)  and  (2.52)  are  given  by  (for  details,  see 
[BSW],[BSW2]) 

«  =  -SiA0)j3NI, 

=  -—dMfl. 
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3  Semigroup  Formulation  and  Well-Posedness 


3.1  System  Formulation 

Let 

V  =  //'(th,fl,)  x  H2o(e1,02)  x  [ll^(0u02)\'. 

and  introduce  an  inner-product  on  V  by 

/  (  Vl\  (  V2\\ 

{  I  )  =  ((6  -  adi)dui,3u2) 

X\zJ  \z2)/v  (3.1) 

+  (&d2wli S2ui2)  +  ((1  —  a)AdZi,  dZ2 ), 

where  (■.■■)  is  an  L2 -inner-product  with  respect  to  Lebesgue  measure.  We  require 

£1  —  aa  1  >0  1  —  a  >  0 

in  (3.1)  in  order  for  this  to  define  an  inner  product  for  V. 

Let 


7i  —  L2(6,i,6I2)  x  L2 ($1 ,  O2 )  x  [L2(6h,  #2)]4. 


Introduce  an  inner-product  on  7 i  by 
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Next,  let 
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We  introduce  an  inner-product  on  LV  by 
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We  regard  the  operators  A  and  B i  given  in  (2.83)  and  (2.84)  as  elements  of  £(  V,  V*). 
Define 
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3.2  Abstract  Cauchy  Problem 


We  can  verify  that 


{Vr],r])w  <  0 


provided  that  g(—r )  >  0  and  ^  >  0. 


Let 


x  =  (v,  w,  Z,  v,  w,  Z,  Tj(t  +  ■, -))T 

■v  =  (V,  |  ■  Iv  X  (7Y,  I  •  Ik)  X  (W,  I  ■  |w). 
Thus,  we  have  the  formal  abstract  Cauchy  problem 


x(t )  =  Ax(t)  +  F  in  X, 


(3.4) 


where 


with 
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A  = 
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v) 

dom  A 


(«•  <M)r  e  A’ 


v/i  £V,i)£  dom  'Pi 
—  /!</>  +  Kr]  £  Ft  J 


Remark  3.1  In  (2.78)  the  underlined  terms  have  coefficients  of  the  same  order  in 
/t/ii  if,  for  example,  r  +  T  +  A  hj\fR.  The  term  dMg/R  is  the  sum  of  shear  stresses 
in  the  elastic,  constraining,  and  piezo  layers.  Since  these  layers  are  considered  elastic, 
the  contribution  of  dM/R  to  the  in-plane  forces  in  these  layers  is  considered  negligible 
and  ignored  (see  [BSW],  [M],  [BSW2]).  Thus,  the  underlined  terms,  of  which  dMsjR 
is  one,  are  ignored  in  what  follows.  In  subsequent  work  we  will  consider  these  terms  as 
perturbations  and  determine  the  extent  of  their  effect  on  the  dynamics  of  the  beam. 


We  now  proceed  with  Remark  3.1  in  mind.  We  observe  that  A  as  defined  above 
is  the  restriction  of  A  £  £(  Jh  3^*),  where  y  =  V  X  V  X  VV,  defined  by  the  sesquil  in  ear 
form:  cr^,  T)  =  (*4.$,  so  that  (t($,  T)  =  {*4.$,  T ) x  for  $  £  dom  A,  'ih  £  X . 

We  can  verify  that 
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for  some  positive  constant  A0.  Here  i/>  can  be  thought  of  as  a  constant  function  in  W. 
Now,  it  is  clear  that  A  —  Ao /  is  dissipative. 

Let  (/,  g,  h)T  £  X  and  consider  the  equation 
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(A  -  A) 

for  A  >  Aq.  This  equation  is  equivalent  to 
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(f 
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(3.5) 


\(f>  —  if;  =  f 

\if)  +  A<f>  —  Kr]  =  g 
—  y'  +  (A  —  T>)r]  =  h. 


(3.6) 


Rearranging  and  solving  for  r/  in  terms  of  other  variables  we  obtain 


AV  +  M  -  (K(  1  -  eAs)</>)  =  M(A  -  V)  '(I,  f)\  ■  g  ■  A/.  (3.7) 

Define  a  sesquilinear  form  ((/>,  ip)  on  V  X  Vby 

=  A2(^,  ip)  +  (Ao,  <n  -  (K{  1  -  eXs)(j),ip). 

It  is  readily  seen  that  g,\  is  coercive.  Using  Lax-Milgram’s  lemma  we  solve  (3.7)  for 
(p.  In  turn  we  have 

'tp  =  \<f>  —  f  £  dom  A  =  V 

t]  =  (A  —  £>)_1(ft  +  0)  (E  dom  'D 
—^4  +  Krj  =  g  —  \ip  G  74. 

Thus,  we  can  solve,  (3.6)  for  (</>,  ^,r/)T  (E  dom  >1.  Thus,  „4  generates  a  Co-semigroup 
/'(/.)  on  X  which  satisfies  |7'(/.)|  <  l/C  . 

The  above  theory  is  adequate  to  treat  the  piezo  embedded  constrained  layer  curved 
beam  of  Section  2  so  long  as  there  are  no  input  voltages  to  the  patches  and  the 
normal  external  forces  qn  are  sufficiently  smooth  (i.e.,  in  i2(ffi,$2))-  However,  when 
the  patches  are  activated,  the  forces  qg  and  dmx  are  not  in  $2)  (in  this  case 

(qg:dmx)  is  in  V*)  and  a  theory  that  extends  the  sense  of  equation  (3.4)  and  the 
semigroup  T{t)  to  a  larger  space  containing  y*  is  necessary.  The  ideas  to  carry 
out  this  extension  are  by  now  rather  straightforward  although  technically  somewhat 
tedious.  The  generic  arguments  for  this  procedure  (often  referred  to  as  “Haraux 
extrapolation”)  are  given  in  [BKW]  for  general  examples  and  specific  arguments  for 
second  order  systems  written  as  first  order  abstract  Cauchy  systems  (such  as  the 
example  of  this  paper)  are  given  in  detail  in  [BSW2],  This  leads  directly  to  an 
extension  of  T(t )  on  X  to  a  Co-semigroup  T{t)  on  a  space  y *  containing  y*  = 

V  X  V*  X  L2[—  r,  0;  V*)  where  we  have  assumed  the  usual  Gelfand  triple  framework 

V  ‘H  'H~  ‘  >  V%  and  y  ^  X  ^  V*  where  y  =  V  X  V  X  L2(-r,  0;  V)  and 
X  =  V  X  7i  X  L2[— r,  0;  V).  Thus,  the  theory  above  is  readily  extended  to  treat  the 
general  case  of  active  patches  in  the  models  of  Section  2.  We  shall  not  give  the  details 
here. 


4  Approximation 

The  equations  (2.78),  (2.79),  and  (2.81)  represent  an  infinite  dimensional  dy¬ 
namical  system,  and  thus  whether  interested  in  control,  parameter  estimation,  or 
just  simply  simulation,  we  desire  a  finite  dimensional  computational  algorithm.  For 


this  purpose,  we  need  an  approximation  theory  to  serve  as  basis  for  computational 
schemes.  In  this  section  we  describe  one  class  of  approximations  that  provides  a  good 
basis  for  computations. 


4.1  The  Approximation  Spaces 


After  considering  the  function  spaces  in  Section  3.1  we  have  chosen  linear  B- 
splines  and  cubic  A-spI ines  to  approximate  the  tangential  and  the  radial  motions, 
v(-,6)  and  w(-,6)  respectively.  We  shall  use  N  to  denote  the  number  of  subintervals 
used  in  partitioning  [^1,^2]-  We  take  the  standard  splines  (see  the  Appendix  of  [BK] 
as  well  as  references  [Sh],  [Sell])  corresponding  to  the  partition  {0,}AO  and  we  obtain 
N  +  1  linear  splines  and  (TV +  3)  cubic  splines.  Looking  at  the  function  spaces  chosen 
we  see  that  we  obtain  N-\-l—2  —  N  —  l  linear  splines,  and  iV  +  3  —  4  =  JV  —  1  cubic 
splines.  We  denote  these  linear  and  cubic  splines,  {CltE  respectively. 

Using  them  as  basis  elements,  we  define  the  following  space 


KN-l  N —1  N  —  1 

E  E  AC  E  C  C 

4=1  4  =  1  4  =  1 


N- 1 

■■■>  E 

4=1 


a- 


2,4 


«4,  af,  Pi  £  H 


(4.1) 


For  each  N ,  V  x  C  V  C  H.  and  the  sequence  {V'C  provides  an  approximation  space 
to  V.  Letting  .s(v/  =  — jr/M ,  j  =0,1,...,  M ,  set 


Let 


sf  <  s  <  sfij 
otherwise 


(4.2) 


(  M  N- 1  M  N- 1 

wM,N  =  (EE  PyC,(%f  (*),  E  E 

\  j=l  i= 1  j= 1  i= 1 

M  TV  — 1 

J2  £  <E  ft) 


(4-3) 


J  =  1  4  =  1 


As  a  result,  the  approximation  spaces  to  the  state  space  X  are  given  by 


XM’N  =  VN  X  VN  X  WM'N. 


Let 


P%  :H^V 


N 


Pv  =V 


V 


AT 


V, 


M.M 


w  w 


denote  the  respective  orthogonal  projections.  Thus,  for  x  =  £  X,  we  have 

xM,N  =  (<f)N ,  y',A ,  rjM,N)  £  XM,N,  with  pN  =  Py  </>,  tpN  =  P-^p,  and  rjM,N  =  Pq’^tj. 
In  particular 


vM'N  = 


e  n 


M  N- 1 

v?’N  =  e  E  (%f  (*)> 

i= 1  *=1 
M  N- 1 

=  E£#E%?(4  *&€* 

|=i  «=i 
M  N—l 

# N  =  EE4^%f(^  4-^4. 

j=i  *=i 


We  define  the  operator  DM,N  :  W 


M,N  .  .  -\a-M.N 


M.NM.N 


DM’"ri 


*  WM’n  by 

M  /  ! 


E  E  v  «,-!  -  itj)  x"M&uW. 


*'=i  j=i 


iv-l  M  (i# 

EE7*fWW(4,-i). 


*=i  j=| 


iV-l  M  M  \ 

E  E7(4,-43*?MW 

*=1  7  =  1  / 


with  r/E  =  t]2i  0  =  0,  r/f)0  =  (0,  0,  0,  0)T. 

We  can  now  define  the  approximate  system  dynamics  state  operator  A  W’  '  on 

dom  AM’N  =  {{4F,*I>N,tiM’*#'  e  XM’N  :  r,M’N  G  W^} 


^M,iV  +  iTryWAf 

^ M,N  J  ^  + 


4.2  Convergence  Analysis 

Lemma  4.1  The  operator  T>M,A  :  WM,A  — +  >VM,JV  is  dissipative. 


Proof.  Let 


JV-l 

sf  =  E  *  =  1.2 

4  =  1 


4  =  E«W 


«*%"’")  =  {«*  E"!1  E"i  -  rtjKiX?, 

a‘Efe'Ejli^Xi*f> 

We  thus  have 

I°-r9^)(dkDM^p^N,dkp^N)dS  <  1  (kf  1 1  /  .a  -  kill)  f-r9(s)*;?(*)d*  <  0. 


We  also  note  that 


/  g(s){dDM’NAr,^N,dr,^N)ds 

J —r 

=  /  g(s)(dDM,NVArj^[,N ,dVAr)%I,N)ds  <  0. 

J —r 


Thus,  DM’N  is  dissipative. 
Let 


LV 


M,N  _ 


where 


M  N- 1  M  N- 1 

E  E^E*)iTt4E  E  ^b?ti(6)Ef(s% 

"  j= l  «=i  j= l  ?— l 

M  JV-l  x  ^ 

E  E  ■■■■  9ijbi,ii6))  E,('s) )  N’ Si}* ry  €  n  k 

j=i  *=i  y  J 


_M 


Ef{s)  =  < 


Define  an  isomorphism 
by  the  formula 


[  M(s-  sfl^/r, 

*i+l  <  *  < 

i 

rH 

Si-i 

03 

5 

sf1  <  s  <  , 

l  0, 

otherwise 

im,n  ■  yy  —>■ 

WM,N 

M  JV-l 

M  JV-l 

M 


M  N- 1 

*»,»  I  £  £  £  £  (44.. '  • '  •  44)  JJ 

u=i  «=i  i=i  »=i  j=i  *=i 


M  JV-l 


M  JV-l 


M  JV-l 


=  £  £  kHW.Z  E  M‘Kxf . £  £  (4* 44 


.JV  f  M 

Aj  j 


v.;  -  t  ;  i 


j=i  *=i 


J  =  1  S  =  1 


p<,o  =  gji0  =  Dp  =  0.  For  ihM’JV  <E  >VA/A',  DuiM’iV  =  DM’NiM,N^M’N- 
Lemma  4.2  For  A  >  0,  (XI  -  DM’N)~1P^’N  h  -  (A/  -  D)_1h,  h  €  >V. 

Proof.  The  proof  of  this  lemma  proceeds  exactly  in  the  same  manner  as  in  [BMZ],  We 
refer  the  reader  to  that  paper  for  details. 

The  convergence  scheme  presented  is  adequate  to  treat  approximate  solutions  of  the 
homogeneous  version  of  (2.78),  (2.79)  and  (2.81).  To  obtain  convergence  for  solutions  when 
the  right  hand  side  of  this  system  is  in  T* ,  one  must  extend  the  convergence  scheme  to  the 
extended  semigroup  T(t)  on  y*  that  was  discussed  at  the  end  of  Section  3.  The  necessary 
convergence  scheme  will  be  presented  elsewhere. 


5  Concluding  Remarks 

We  have  presented  a  general  functional  analytic  framework  that  can  be  used  to  treat 
a  wide  class  of  curved  active  constrained  layer  structure  models.  This  work  extends  our 
previous  work  ([BMZ])  where  we  did  not  consider  the  shear  contribution  in  the  viscoelastic 


layers  and  constraining  layers  were  absent.  The  inclusion  of  shear  effect  and  the  constrain¬ 
ing  layers  makes  the  model  more  realistic.  We  have  illustrated  how  to  treat  one  type  of 
viscoelastic  hysteresis  in  these  models.  Moreover,  an  approximation  technique  that  can  be 
used  as  a  foundation  for  computational  methods  has  been  presented. 

The  ideas  presented  here  are  quite  general  and  can  be  used  (with  minor  modifications 
and/or  extensions)  to  treat  most  of  the  active  constrained  layer  structure  models  found  in 
engineering  and  scientific  research  literature. 
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first  prompted  our  interest  in  ACL  curved  structures  such  as  that  modeled  in  this  and  our 
previous  paper.  They  are  also  indebted  to  Dr.  Beth  Munoz  and  Dr.  Lynn  Yanyo  of  the 
Thomas  Lord  Research  Center,  The  Lord  Corporation,  for  fruitful  discussions  related  to 
viscoelastic  materials. 


6  Appendix 


We  present  here  the  matrices  Ai,  B  i  used  in  (2.83).  In  (2.78),  (2.79),  Mg,  the  sum  of  the 
internal  moments  of  the  elastic  constraining,  and  piezo  layers  about  the  neutral  line  of  the 
elastic  core  is  approximated  by  —  'd2  o  £2  0  d2io,  where  £2  is  as  in  (2.46)  and  it  is  understood 
that  derivatives  are  taken  in  the  distributional  sense  where  necessary.  Some  of  the  entries  of 
the  matrices  need  explanation,  and  we  do  that  immediately  after  presenting  the  matrices. 
They  are  given  by 

/  ^11  zfi3\ 

Ai  =  I  .I21  A22  A23  I 

V  A32  A33  / 

where 


An  =  -ldohod  +  libR-'a 2  +  (A,  Ji)\ 

An  =  -bJj^d  o  6  -  HR  '(n-2  +  (hj,  J2))d  +  5  (d(0),  J2)  d2 -^do&od2-  ±rlpd  °£2°d2, 
Ais  =  !<£>!,/■)-  ld(D{0),-), 

An  =  ^(bR-H  2  +  6)3, 

A22  =  i (bR-1^  +  D282)  +  lpbR-H2  o  ^  o  d 2, 

A23  =  ±(D(0)-D1,.), 

-431  = 

zf  32  =  +  ^77-^3; 

A33  =  ~^d2  - 
and 


Bi  = 


B  n 

-012 

B 

-021 

022 

0 ; 

-031 

032 

0 

where 


Bn  =  j(aid2  -  a2), 

Bn  =  '-bit  >2  -  ai)d  -  UR-1  pe,  j2)  d 2  +  Ur  'm\ 

Bn  =  UR-1  (J^,d-) +  1(^,1-), 

B2i  =  jiai  +  02)9, 

B22  =  -UR-U^^-D2, 

B23  =  j  {D\  —  bR~1Jve,  ■)  , 

B31  =  -pE4, 

B32  =  ^£"5  +  ^-£-6, 

£^33  =  ^92  o  A  +  ^T-P. 

We  now  explain  the  entries  in  A\,B4  as  needed.  From  (2.82)  recall  dV(t,6)  =  Z(t,8). 

In  matrix  ,lp.;  {l)\.  /•)  ( Z)  -  where  (•,•)  stands  for  inner  product  in  7v4.  We 

remark  that  A\  and  B i  act  on  the  vector  ( v ,  u),  zv)T.  Continuing, 

W)  =  {D(6)-DhZ) 

p!  -  bR~1Jve,-)(Z)  =  p!  -  bRT1Jve,z'S} 

H:n  =  l-{i-Vve){dv)PAT-1j1 

E2(w )  =  ( ()2  u- 1  h 

i  =  -1  pe(d2w)J2. 

Next,  we  deal  with  matrix  /F| .  In  matrix  B\, 

(Jve,  d-)  (Z) 

(Di,E)(Z) 

{/;,  -  bit  \r.,-)(Z\ 

E4(v) 

E5(w ) 

£^e(w) 
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